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Abstract 

In this paper we describe General Covariance Union (GCU) and show that 
solutions to GCU and the Minimum Enclosing Ellipsoid (MEE) problems are 
equivalent. This is a surprising result because GCU is defined over positive 
semidefinite (PSD) matrices with statistical interpretations while MEE involves 
PSD matrices with geometric interpretations. Their equivalence establishes an 
intersection between the seemingly disparate methodologies of covariance-based 
(e.g., Kalman) filtering and bounded region approaches to data fusion. 



1 Introduction 

Positive semidefinite matrices are often chosen to represent uncertainty in state es- 
timation and control algorithms because of their special linear-algebraic properties. 
The Kalman filter, for example, uses PSD matrices to reprepresent covariance upper 
bounds on the second central moments of unknown probability distributions relating 
to the state of a system of interest. Bounded region filters, by contrast, use PSD 
matrices to define ellipsoidal regions which bound the state of the system of interest. 

The choice between representing uncertainty with covariance upper bounds versus 
ellipsoidal bounded regions leads to very different data fusion algorithms with very 
different filtering and control properties. The statistical interpretation associated 
with covariance matrices leads to use of the Kalman fusion equations in the case of 
independent estimates or Covariance Intersection (CI) when the statistical relation- 
ships among estimates to be fused cannot be established [9]. Under their respective 
assumptions, Kalman and CI yield a fused estimate with the smallest possible co- 
variance that is guaranteed to be an upper bound on the second central moment of 
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the unknown probability distribution defining the state of the system. A bounded 
region filter, on the other hand, determines the minimum-size elhpsoid that bounds 
the intersection of the given elhpsoids which are assumed to bound the true state of 
the system. In summary, one framework is statistical while the other is geometric. 

Despite the use of PSD matrices and the exploitation of linearity properties, covari- 
ance and ellipsoidal frameworks employ very different mathematical techniques and 
have almost completely disjoint literatures. (This relationship also holds true more 
generally between statistics and computational geometry.) A consequence of this fact 
is that the sophisticated tools developed separately within each framework seem to 
have no applicability within the other framework. In this paper we make progress 
toward bridging this gap by establishing a surprising equivalence result. 

The structure of this paper is as follows: The next section provides background for 
the Kalman filter, Covariance Intersection (CI), and Covariance Union (CU). We then 
present General Covariance Union (GCU) and show that its solution is equivalent to 
determining the Minimum Enclosing Ellipsoid (MEE) of the ellipsoids defined by its 
covariance matrix arguments. In other words, GCU and MEE are equivalent. We 
end with a discussion of the implications of this unexpected result. 

2 The Covariance Upper Bound Framework for 
Data Fusion 

An estimate of the squared error, or covariance, associated with measurements from 
a particular sensor can be modeled empirically by examining measurements taken 
of reference objects whose true states are known. This permits an error covariance 
matrix to be associated with subsequent measurements of objects whose true states 
are not known. For example, the measured position of an object in two dimensions 
can be represented as a vector a consisting of the object's estimated mean position, 
e.g., a = [x,y]^, and an error covariance matrix A that expresses the uncertainty 
associated with the estimated mean. If the error in the estimated mean vector is 
denoted as a, then the error covariance matrix is an estimate of the expected squared 
error, E[aa^]. 

Ideally, A would equal E[aa^], but this is not generally possible to achieve in practice 
because measurement and process models are never perfect. To accommodate the 
effect of model errors, prediction and measurement covariances are typically overes- 
timated so as to avoid underestimating the actual squared errors. In other words, a 
more conservative overestimate of errors is deemed preferable to underestimating the 
errors. One of many reasons for this preference is that it avoids the consequences of 
spuriously small errors causing a covariance matrix to become singular or numerically 
unstable. 



2 



Formally, a mean and covariance estimate is said to be consistent (or conservative) if 
and only if A - E[aa"^] is positive semidefinite (i.e., has no negative eigenvalues): 



A y Efaa 



The statistical properties associated with a mean vector and its associated covariance 
upper bound are completely defined by the definion of consistency: all that can be 

said about a mean and covariance pair, (a, A), is that the covariance matrix. A), is 
greater than the expected squared error in its associated mean, a). 

Having established a rigorous definition of what constitutes a consistent or conser- 
vative estimate, it is possible to certify the performance of the Kalman filter and 
Covariance Intersection. 



2.1 The Kalman Filter 

Given two mean and covariance estimates (a. A) and (b, B), the data fusion problem 
of interest in this paper consists of determining a fused estimate (c, C) that is guar- 
anteed to be consistent and summarizes the information in the two estimates with 
error (in terms of the size of C) that is less than or equal to that of either estimate. If 
the two estimates are consistent and presumed to be statistically independent, then 
a joint estimate can be constructed as: 



Letting a and b denote the errors in the respective mean estimates, the key property 
of the joint covariance estimate is that it satisfies: 

" A 
B 

where the RHS matrix represents the true but unknown joint error covariance, which 
has zero cross covariance, ab-^ = 0, due to the assumption of statistical independence. 
The estimated joint covariance is a conservative estimate of the true joint covariance 
because in practice A >z E[aa"^] and B y E[bb"^]. The latter inequalities hold by 
design in that intentional efforts are made to ensure that estimate error covariances 
do not underestimate the actual squared errors associated with sensor and kinematic 
models. 

Given a consistent joint covariance for two given n-dimensional estimates, the Kalman 
filter defines the optimal linear projection of the 2?7,-dimensional joint estimate back 
to the n-dimensional state space of interest. The result of the Kalman projection is 
a mean and covariance estimate (c, C) that represents the optimal fusion of the two 
given mean and covariance estimates. In fact, if there is no additional information 





A 




B 



(2) 



E[aa^i 






E[bb^ 



(3) 
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available (e.g., distribution information), then the Kalman fusion estimate is optimal 
according to virtually any error criteria |7]. 

In the case of statistically independent estimates (ai, Ai), (a2, A2), (a^, A^), the 
Kalman fusion equations have a particularly simple form [7]: 

C = {A-' + A,' + ... + A;,y 

c = C(Ar^ai + A2^a2 + ... + A„^a„). (4) 

If its underlying assumptions hold (i.e., consistency and independence), then the 
above Kalman equations ensure that the fused estimate (c, C) is consistent, and C ^ 
Aj Vi, 1 < z < m. However, any presumption of statistical independence in practical 
data fusion contexts should be carefully considered. Specifically, virtually any sensor 
is subject to time-correlated errors induced by the particular conditions of its use 
(e.g., changes in temperature, platform vibrations, relative humidity), and errors 
associated with the nonlinear transformation of its measurements (e.g., from local 
spherical coordinates to a global coordinate frame) are deterministic and therefore 
non-independent. 

If estimates (ai, Ai), (a2, A2), (a^, A^) are each consistent, but not completely 
independent, then it is possible for the Kalman fused estimate (c, C) to be inconsis- 
tent. In fact, if Aj = a^af Vi, 1 < z < m, then any degree of correlation guarantees 
inconsistency, i.e., C ^ E[cc'^]. The key point is that a Kalman fused estimate is not 
guaranteed to be consistent even if each of its given estimates are consistent. The 
reason why the Kalman filter fails is because although two given estimates (a, A) 
and (b, B) may be individually consistent, the implicit joint covariance may not 
be if independence is assumed when the cross covariance between the estimates is 
X = E[ab^] ^ 0. Specifically: 

" E[aa^] 
E[bb^] 

In other words, the Kalman filter only fails to produce a consistent fused estimate 
when the implicit jomt estimate is inconsistent. Although the equations are very sim- 
ple and elegant for independent estimates, the Kalman filter is also defined generally 
for any consistent joint covariance with X 7^ [7]. Therefore, the Kalman filter is 
guaranteed to produce consistent estimates as long as the given estimates are con- 
sistent and their cross covariance is known [3]. Unfortunately, this poses significant 
challenges. The first challenge is that the cross covariance information must in prin- 
ciple be determined exactly. This can be seen by examining the difference between 
two joint covariance matrices with different cross terms: 

A X 
X^ B 

The difference matrix is not PSD for any case in which X 7^ Y. The need for 
absolutely perfect cross covariance information presents difficulties when estimates are 



E[aa^] X 
X^ Efbb^l 



(5) 
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Y 








B 





X- Y 

(X - Y)'^ 



(6) 



4 



the products of nonlinear operations (e.g., coordinate transformations, kinematic time 
projections, human-derived estimates) because the error processes are not perfectly 
modeled. For example, the same approximate nonlinear transformation equations 
may be applied to convert different radar observations of an object to a common 
coordinate frame, so the errors committed are clearly not independent, but it may 
not be possible to determine exact cross covariances. 

Kalman's original derivation of his eponymous filter was based on orthogonal pro- 
jection theory, and the fact that there exists a simple Bayesian interpretation of the 
result when error distributions are Gaussian was presented only as an interesting 
special case 0. However, many subsequent references to the Kalman filter incor- 
rectly suggest that the Kalman filter requires assumptions of Gaussianit}0. It turns 
out that the assumption of estimate independence is actually the only problematic 
assumption because it typically cannot be guaranteed in practice. Relaxing the in- 
dependence assumption leads to the more general fusion equations of the Covariance 
Intersection (CI) method which has proven invaluable for a wide variety of practical 
apphcations[2l H]. 



3 Covariance Intersection (CI) 

The general mean and covariance data fusion problem can be formulated in terms of 
the joint covariance structure that implicitly exists between a given pair of estimates 
(a. A) and (b, B): 

A X 
B 

where X represents the actual, but unknown, cross covariance between the two esti- 
mates. If X were known, then it would be possible to apply more general formulations 
of the Kalman filter equations to produce an optimal fused estimate. Unfortunately, 
these generalizations only guarantee consistency if the cross covariance is known ex- 
actly, i.e., it cannot be conservatively approximated in any way analogous to the way 
conservative covariance estimates are used. 

Without knowledge of X, the only way to ensure consistency in the application of the 
Kalman filter is to identify a joint covariance that is guaranteed to be consistent based 
on the information available. In the present context, therefore, a joint covariance 

^In fact, one commonly-cited motivation for investigating the application of neural networks and 
fuzzy logic is the claim that the Kalman filter imposes restrictive Gaussianity assumptions that 
cannot be satisfied in many applications. The fact is that the use of covariance upper bounds was 
recognized as necessary when the first Kalman filters were implemented in the late 1960s, and such 
bounds are incompatible with PDF interpretations. It was shown by Jazwinski in his classic 1970 
book that the standard practice of using covariance upper estimates does not undermine the integrity 
of the Kalman filter [3]. 



(7) 
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matrix M must be determined such that: 



M ^ 



A X 

X^ B 



for every possible cross covariance X. It can be inferred from the symmetry of the 
unknown cross covariance information (i.e., M must be consistent for any instanti- 
ation X = A' and for X = —X) that the off-diagonal blocks of M should be zero, 
and its diagonal blocks must be sufficiently larger than A and B to account for the 
effects of all possible degrees of correlation among the error components of the mean 
estimates a and b. 

It has been shown (appendix 14 of [9]) that a consistent and tight joint covariance 
M can be generated by selecting a scalar value u, < u < 1 as: 



iA 








B 



A X 
X^ B 



(9) 



where u is chosen to minimize the size (e.g., determinant) of the covariance produced 
by the Kalman filter update equations for the estimates (a, ^A) and (b, ^z^B). This 
can be generalized for an arbitrary number of estimates: 













-A, 












Ai Xi,2 . 

-^2,2 A2 
Xm 1 Xm 2 



X2,m 



(10) 



where Yl^i = 1) the parameters are determined to minimize the covariance 
resulting from the fusion of the n estimates. The above inequality provides a general 
(and optimal) mechanism for obtaining a consistent joint covariance when given only 
the block diagonals of an unknown joint covariance matrix. The use of appropriate 
w-paramterized covariances in the Kalman update equations yields the CI fusion 
equations. 

The general approach of determining consistent joint covariances can also be applied 
to solve related problems. For example, given estimates (a, A) and (b.B) that are 
correlated to an unknown extent, the covariance of a -|- b can be computed as ^A + 
Y^B. This is referred to as Covariance Addition (CA). 



3.1 Covariance Union (CU) 

Covariance Intersection addresses the general form of the data fusion problem for 
mean and covariance estimates, but in practice a different problem can arise before 
data fusion can even be performed. Specifically, what is to be done if two estimates 
(a. A) and (b, B), purportedly relating to the state of the same real-world object, are 
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determined to be mutually inconsistent with each other, i.e., the differences between 
their means is much larger than what can be expected based on their respective error 
covariance estimates? For example, if two mean position estimates differ by more 
than a kilometer, but their respective covariances suggest that each mean is accurate 
to within a meter, then clearly something is wrong. 

One mechanism for detecting statistically significant deviations between estimates is 
to compute Mahalanobis distances [6]. The Mahalanobis distance between estimates 
(a. A) and (b, B) is defined as: 



which is essentially just the squared distance between the means as normalized by 
the sum of their respective covariances. Intuitively, if the covariances are large, then 
a large difference between the mean vectors a and b is not surprising, so the Maha- 
lanobis distance is small. However, if the covariances are very small, then even small 
differences between the means may yield a large Mahalanobis distance. A large Maha- 
lanobis distance may tend to indicate that the estimates are not consistent with each 
other, but a user-defined threshold is required to define what constitutes an accept- 
able deviatioiJl. When the threshold is exceeded, the estimates are regarded as being 
contradictory and some kind of action must be taken. Resolving such inconsistencies 
among estimates is sometimes referred to as deconfliction 



The Covariance Intersection method guarantees consistency as long as the estimates 
to be fused are each consistent. In the deconfliction problem it is only known that 
one of the estimates, either (a. A) or (b,B), is a consistent estimate of the state of 
the object of interest. Because it is not generally possible to know which estimate 
is spurious, the only way to rigorously combine the estimates is to form a unioned 
estimate, (u, U), that is guaranteed to be consistent with respect to both of the two 
estimates. Such a unioned estimate can be constructed[10j by computing a mean 
vector u and covariance matrix U such that: 



where some measure of the size of U, e.g., determinant, is minimized. This Covariance 
Union (CU) of the two estimates can be subsequently fused with other consistent 
estimates using CI. The CU equations generalize directly for the case of m > 2 two 

^It must be emphasized that the use of a threshold on Mahalanobis distance is not the only 
possible mechanism for identifying potentially spurious estimates, but some user-defined mechanism 
is required. Otherwise there is no way to distinguish fault conditions from low probability events. 
In other words, models for fault conditions are inherently application-specific. 
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estimates: 



U h Ai + (u-ai)(u 
U h A2 + (u-a2)(u 




U h A^ + (u-a„)(u 




(13) 



Intuitively, the CU equations simply say that if the estimate (a, A) is consistent, then 
the translation of the vector a to u will require its covariance to be enlarged by the 
addition of a matrix at least as large as the outer product of (u — a) in order to be 
consistent. The same reasoning apphes if the estimate (b, B) is consistent. Covariance 
Union therefore determines the mean vector u having the smallest covariance U that 
is large enough to guarantee consistency regardless of which of the two given estimates 
is consistent. The resulting covariance may be significantly larger than either of the 
given covariances, but this is an accurate reflection of the actual uncertainty that 
exists due to the conflict between the two estimates. The key fact is that the CU 
estimate satisfles the deflnition of consistency. 

As a simple example of a CU construction, consider two estimates (a. A) and (b, B) 
of the location of an object observed from two nodes in a network. The estimate from 
the first node places the mean position at a = [0, 0]^, and the second node places it 
at b = [4, 4]"^, and each has an error covariance equal to the identity matrix I. If it is 
determined that the two estimates are statistically inconsistent with each other, thus 
implying that one of the estimates is not a consistent estimate of the object's location, 
then deconfliction must be performed. The optimal CU deconflicted estimate is: 



It is straightforward to verify that this estimate (u, U) is in fact consistent with 
respect to either/both of the estimates (a. A) and (b,B). If (a. A) is a consistent 
estimate of the target's state, then the covariance U for mean u must be greater than 
or equal to A + (u — a)(u — a)^, which it is. It can be verified that the estimate 
(u,U) is similarly consistent with respect to the estimate (b,B). Therefore, if either 
of the two estimates represents a consistent estimate of the state of the object, then 
the CU estimate is also consistent. 



Consistency of the CU equations rests on an implicit assumption that the estimates 
to be combined are not correlated. Specifically, the CU inequahties for the estimate 



4 General Covariance Union (GCU) 
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(u,U): 



U y Ai + (u - ai)(u - ai)^ 
U >z A2 + (u-a2)(u-a2)^ 

U y A„ + (u-a™)(u-a^)^ (15) 

technically hold only if the errors associated with Aj are uncorrelated with a,. This is 
true for most types of common process models, but not for certain kinds of recursive 
control and colored noise models. 

For example, suppose that we are given a ID estimate (0, 0), i.e., zero mean and zero 
variance. The underpinning assumptions of CU imply that if the mean is translated 
to 1, the new consistent estimate must have variance + 1^ = 1. If instead the mean is 
translated to 2, the variance must be 0+2^ = 4. However, suppose a 2-step translation 
is performed: First the estimate is translated 1 unit to produce (0 + 1, + 1^) = (1, 1), 
then that estimate is translated 1 more unit. According to CU the result should be 
(1 + 1,1 + 1^) = (2, 2), but this is clearly incorrect because translating the estimate 
(0, 0) by 2 units must produce (2, + 2^) = (2, 4) to ensure consistency. The problem, 
of course, is that the translations in the sequence of steps are not independent. 

More generally, given a mean and covariance estimate (a, A), the covariance of a+x is 
equal to A + xx^ if and only if the estimate error, a, is independent of the translation 
X. If they are correlated to an unknown extent, then the covariance of a + x must be 
formulated using Covariance Addition (CA) as ^A + j^xx^. Applying CA to the 
CU inequalities (ITS!) gives the generalized GCU formulation: 

u y 
u y 



U y — A„ + - — - — (u - a„)(u - a^)^ (16) 

where the optimization problem now involves n free parameters (0 < cuj < 1), like CI, 
but with the semidefinite inequalities of standard CU. This appears to represent a 
daunting challenge to solve efficiently. However, plots of solutions for low- dimensional 
problems provide a critical insight. Specifically, Figure 1 depicts the 1-sigma solu- 
tion contour that minimally encloses the 1-sigma contours of the estimates that are 
unioned. This suggests - but of course does not prove - that the GCU solution is 
equivalent to the corresponding MEE solution obtained by interpreting the covari- 
ances as ellipsoids defined by their 1-sigma contours. That they are in fact equivalent 
is established in the following section. 



— Ai + (u - ai)(u - ai)'^ 

UJi 1 — Ui 

—As + — ^ — (u - a2)(u - as)^ 

U2 I - UJ2 
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Figure 1: A typical General Covariance Union (GCU) optimization problem. The five small covari- 
ance ellipses represent the five estimates to be unioned. Surprisingly, they seem to be circumscribed 
by the minimal-determinant GCU solution. 



5 The Equivalence of GCU and MEE 

This section provides a formal proof of the equivalence of MEE and Our treat- 

ment of ellipsoid enclosure follows the discussion in §3.7.1 of [IJ, so several of its 
equations are repeated here for easy reference. 

In pp, the enclosing and enclosed ellipsoids are represented via parameter triples 
{Ai, bi, Ci} where index refers to the enclosing ellipsoid: 

{x e I x^AiX + 2x^bi + Ci < 0}, i = 0...m, (17) 



AT yO and bjAT^k - a > 



The triples {Ai,bi,Ci} are homogenous so the authors normalize the enclosing ellip- 
soid's parameters via Cq = feg ^feo ~ 1 and then use the S-lemma to represent the 
ellipsoid-enclosure constraint matrix inequality: 



^0 ^0 

b^ b^A^\-l 



Ti 



Ai h 
bT 



Ci 



^0, Ti> 0, 



1 . . . m. 



:i8) 



That inequality is nonlinear so they use a Schur complement argument to expand it: 

^0, Ti>0, i = l...m. (19) 





ho 







-1 


bl 





bo 


-A 







Ti 



A, bi 

hj Q 0^ 





•^To make the exposition concise our proof assumes that the covariance matrices are nonsingular, 
but the equivalence holds generally. 
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The minimum-determinant ellipsoid enclosure can then be posed as a Maxdet [TT] 
optimization: 



mmimize 
subject to 



log det ^ 

^0 



>- 




>o, 


i = 1 


. . .m. 






bo 







' Ai hi 


" 




-1 




- n 


bj Q 








bo 


-Ao 











(20) 



^ 0, 



1. 



, m. 



5.1 Covariance ellipsoid enclosure as a Maxdet 



In this section, we translate the ellipse-enclosure formulation f l20|) from the original 
{A, b, c} triples to covariance/mean estimation parameters {u, U). 

The covariance ellipsoid S{u, U) for an estimate at mean u with covariance matrix U 
is defined as: 

S{u,U) = {xeR''\{x-ufU-\x-u)<l} 

= { X G i?" I x'^U-\ + 2x'^{-U-\) + {u^U-\ - 1) < 0} (21) 

Comparison of f|T7|) and fl2T|) reveals the following correspondences: 

A^U'\ h^-U-^u, c^u^U-^u-l (22) 

Apply fl22|) to the nonlinear matrix inequality f lTSj) - we will use the result in §5.21 
The result is a GCU equation so A^ now denotes a constraint covariance, and its 
mean: 



A: 



-1 



-A- ^a,: 



-Ar^ttif ajA-^a^-1 



^0, Ti> 0, 



Returning to the covariance ellipsoid Maxdet, apply fl22|) to fl20|) : 



minimize log det [/ 
subject to [/ :^ 0, Tj > 0, 







-^/-^ti 
-1 i-U-'uf 



A: 



1 . . .m, 



i = 1 
(23) 

(24) 







-Ar^a.f ajA-'ai-l 




m. 



^ 0, 



1 . . . m. 



The matrix inequalities in are nonlinear, due to the (quadratic) U~^u terms. 
Linearize them with change of variables v = —U~^u. Also, let W = to obtain a 
more standardized Maxdet formulation: 

minimize log det 

subject to W >- 0, Tj > 0, i = 1 . . .m, 



W V 

T 1 T 

V —1 V 

V -W 



- Ti 



:-A-'a,f ajA-'a^-l 




^0, 



= 1 . . , 

(25) 



m. 
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After solving fl25l) for v and W, the optimal {u, U) are recovered as U = W ^ and 
u = —Uv. 



5.2 Enclosure of a single ellipsoid 

Consider a GCU constraint from (fT6l) for a single estimate (a, A) where A is full rank: 

U h -A+ ^—(u-a)(u-af, ue[0,l] (26) 
w 1 — u 



We shall prove that fl2B]) defines the set of all estimates (-u, whose covariance 
ellipsoid S{u, U) encloses S{a, A). It then follows that the combined GCU constraints 
( IT6|) define the set of all covariance ellipsoids which enclose all the £{ak,Ak) so the 
solution to the f l25|) Maxdet formulation is a minimum-determinant solution to (fT6|) . 

Note that fl2^ only depends on the difference between u and a rather than their 
absolute values, and must hold for any arbitrary value of u. Therefore we may, without 
loss of generality, assume that a = (simple coordinate shift) which simplifies ( 126|) 
to: 

^ ue [0, 1] (27) 



u y -A 

UJ 



uu 



l-UJ 

A similar argument can be applied to a single instance from the ellipsoid-enclosure 
inequalities (123|1 : it must hold for arbitrary values of u, and ellipsoid enclosure is a 
geometric property unaffected by coordinate shifts. So we may again assume that 
a = 0, which simplifies ( 123|) to: 



u 



u 



u'U' 



u 





' A-^ " 


— T 


_ 0^ -1 _ 



^ 0, r > 



(28) 



Equation (!28|) implies the following scalar inequality for the lower right-hand main 
diagonal entry: u^U~^u — 1 + r < 0. Since u^U~^u > we must have r < 1. Make 
that restriction explicit, to bring r in line with the equation ( 127|) u parameter: 



Y u^U-^u-l 





' A~^ " 


— T 


_ 0^ -I 



^ 0, TG [0,1] 



(29) 



We will prove equivalence by demonstrating that f l7r|) and f l2I?]) are mutual inverses. 
But first we must consider the singular points. For (!27|) they are at a; = and uj = 1. 



As u approaches 0, U becomes unbounded and therefore £{u, U) encloses any (bounded) 
covariance ellipsoid £(0,^4). The behavior for u = 1 depends on the value of u: if 
M = then the ellipsoids are concentric and (l27|l reduces to the expected enclosure 
requirement U A. If w 7^ then the right-hand side approaches the sum of A 
and an unbounded rank-1 adjustment in the direction of u. The result is an elliptic 
cylinder - with radial axis through and u - which tightly encloses £^(0, ^4). So both 
singular points lead to enclosure of £"(0, A). 
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Now that the singular points have been considered, we can simply invert 0271) : move 
the rank-1 term to the left-hand side and apply the Sherman- Morrison formula. The 
result simplifies to: 

. , {u-'u){u-'ur 



u- 



(30) 



1 — uj — u^U 

The rank-1 term in fl30|) is bounded and positive semidefinite, since it is related to 
the bounded and positive semidefinite rank-1 term ''^w^ from f l27|) . Therefore, 
1 — OJ — u^U^^u > so we can apply a Schur complement: 



1 — OJ — u^U^^u 



Finally, rearrange terms: 





' A-^ " 


— OJ 


_ 0^ -I 



h (31) 



^ 0, we (0,1) (32) 



Equation (!32|) has the same form as (!29l) (except for the singular points a; = and 
UJ = 1 which were considered earlier). This completes the proof. 



6 Discussion 



In this paper we have established a connection between covariance-based and bounded- 
region models of uncertainty. In particular, we have shown that GCU and MEE have 
equivalent solutions. Thus, techniques for solving one problem can be directly applied 
to solve instances of the other. This equivalence is also suggestive of a potentially more 
general framework that subsumes the covariance and bounded-region approaches. 
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